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PACS 47 . 61 . -k - Micro- and nano-scale flow phenomena 

PACS 83 . 10 . Rs - Computational simulations of molecular and particle dynamics 
PACS 87. 15. -V - Biopolymers: structure and physical properties 

Abstract. - The non-equilibrium structural and dynamical properties of a semiflexible polymer 
confined in a cylindrical microchannel and exposed to a Poiseuille flow is studied by mesoscale hy- 
drodynamic simulations. For a polymer with a length half of its persistence length, large variations 
in orientation and conformations are found as a function of radial distance and flow strength. In 
particular, the polymer exhibits U-shaped conformations near the channel center. Hydrodynamic 
interactions lead to strong cross-streamline migration. Outward migration is governed by the 
polymer orientation and the corresponding anisotropy in its diffusivity. Strong tumbling motion 
is observed, with a tumbling time which exhibits the same dependence on Peclet number as a 
polymer in shear flow. 



Introduction. — The flow properties of soft matter 
jilnder confinement liave been studied intensively in re- 
cent years, stimulated by the desire to unravel the be- 
havior of these materials in nano- and microfluidic de- 
vices. In such systems, confinement effects, a spatially 
Varying shear rate, and hydrodynamic interactions due to 
the presence of impenetrable walls, lead to phenomena not 
observed in bulk systems. Often these three effects are mu- 
tually competing, and emergent properties are a result of 
their combined action. 

■ Flexible polymers exhibit interesting conformational, 
dynamical, and flow properties in the presence of walls 
and in channels [l]-[6]. In particular, cross-streamline mi- 
gration away from the wall has been observed pilSHT^. 
which crucially depends on hydrodynamic lift force ex- 
erted by the wall [5J|ni[T3]- In addition, in channel flow 
flexible polymers migrate away from the channel center, 
resulting in a decrease in concentration at the centerline 
[HHJH]. This originates from the large flow-induced poly- 
mer stretching and alignment in the high shear-rate re- 
gion in the vicinity of a surface, which implies a smaller 
lateral diffusion than that of the coiled conformations in 
the channel center. Similarly, migration has been found 
for rods in Poiseuille flow, which has been attributed to 
their anisotropic mobility (parallel and perpendicular to 
the rod axis) ; migration towards the surface has been ob- 
served when wall hydrodynamic interactions are neglected 

[mn]. 



Semiflexible polymers, such as actin-fllaments or short 
fragments of DNA, are neither rodlike nor are they able 
to undergo large conformational changes. This poses ques- 
tions as: What is the influence of hydrodynamics on the 
polymer distribution? Does stiffness influence migration? 
How are the conformations affected by the flow? Exper- 
iments with actin filaments in microchannels provide ev- 
idence for a strong infiuence of flow on semiflexible poly- 
mers [17]; theory and simulations have demonstrated the 
presence of a hydrodynamic lift force in shear flow [T3] . 

The proper account of hydrodynamic interactions (HI) 
is essential in simulation studies of fluid flows in channels 
as is emphasized by the appearance of cross-streamline 
migration. Recently developed mesoscale simulation tech- 
niques, such as Lattice Boltzmann |lllll8j . Brownian dy- 
namics with a hydrodynamic tensor [T] , and multiparticle- 
collision dynamics [TH[5n] , are well suited to study hydro- 
dynamic flows in microchannels and are able to bridge the 
length- and time-scale gap between the solvent and solute 
degrees of freedom. 

In this letter, we employ multiparticle-collision dynam- 
ics (MFC) to study the properties of a confined semiflexi- 
ble polymer exposed to a microchannel flow. We system- 
atically investigate the influence of the flow rate on the 
polymer structure, conformations, dynamics, and distri- 
bution functions across the channel, and provide new in- 
sight into the polymer migration process. Some character- 
istic polymer conformations in flow are displayed in fig. [TJ 
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Fig. 1: Polymer conformations for the Peclet number Pe — 360. 
(a) Aligned and (b) U-shaped conformations in the center of the 
channel, (c), (d) Transient hairpin-like conformations outside 
the center. The flow direction is from left to right. The channel 
diameter is indicated by the black lines |22) . 



In addition, we discuss the polymer tumbling dynamics, a 
property which has been analyzed previously for flexible 
molecules in both shear and microchannel flows [Sldl] ■ 

Model and simulation method . We adopt a hy- 
brid simulation approach to study the properties of semi- 
flexible polymers in flow, where molecular dynamics sim- 
ulations (MD) for the polymer are combined with MPC 
for the solvent MPC is a particle-based simula- 

tion method and proceeds in two steps. In the streaming 
step, the solvent particles of mass m move ballistically for 
a time h. In the collision step, particles are sorted into the 
cells of a cubic lattice of lattice constant a and their rela- 
tive velocities, with respect to the center-of-mass velocity 
of each cell, are rotated around a random axis by an angle 
a P^[^[^ . The fluid is confined in a cylindrical channel 
with periodic boundary conditions along the channel axis. 
No-slip boundary conditions are imposed on the channel 
walls by the bounce-back rule and virtual wall particles 
[lOlllS], and flow is induced by a gravitational force (mg) 
acting on every fluid particle. 

The linear polymer is comprised of point-like 
monomers of mass M, which are connected by linear 
springs with an equilibrium bond length b. Excluded- 
volume interactions arc taken into account by the trun- 
cated 12 — 6 Lennard- Jones potential, with a character- 
izing the bead size and the energy parameter e [27] ■ The 
bending potential 

N-l 

C/, = |^^(i?,+i-H,f (I) 

1=1 



is added to account for stiffness [1D|, where the R4S are 
bond vectors and the bending rigidity Kf, = LpksT/b^ 
is related to the persistence length Lp for the semiflexi- 
ble polymer. Since we consider pressure-driven flows, no 
gravitational force acts on the polymer. 

The interaction of a polymer with the solvent is real- 
ized by inclusion of its monomers in the MPC collision 
step [28]. Between two MPC steps, several MD steps are 
performed to update the positions and velocities of the 
monomers. Extensive studies of polymer dynamics con- 
firm the validity of this procedure [^[?7H^. 

An advantage of the MPC approach is that HI can eas- 
ily be switched off, without altering the monomer diffu- 
sion significantly [501121] • In this case, denoted Brownian 
MPC, each monomer independently performs a stochastic 
collision with a phantom particle with a momentum taken 
from the Maxwell-Boltzmann distribution with variance 
m (Nc) ksT, where {Nc) is the average number of solvent 
particles per collision cell [501131]. 

We employ the parameters a = 130°, h = O.It, with 
r = ma? /kBT {ks is Boltzmann's constant and T is 
temperature), {Nc) = 10, M = m{Nc),b = a = a, the 
fluid mass density g = (Nc) m/a^, ksT/e = 1, and the 
time step in MD simulation Hmd = 5 x 10^'^r. A polymer 
with A'^ = 14 monomers is placed in a cylindrical channel 
of radius R = 8.5a. With the length Lr = {N-l)a = 13a, 
it does not interact with the wall when its center of mass 
is near the channel center. The persistence length is set 
to Lp = 2Na w 2Lr. Hence, the polymer shows rod- 
like behavior at equilibrium. The channel length is 28a. 
Simulations of the pure solvent system yield velocity pro- 
files which agree with the solution of Stokes' equation for 
this geometry. Averages and probability distributions are 
calculated in the stationary state for various independent 
initial conditions. To maintain a constant temperature, 
the velocities arc scaled in every collision step and inde- 
pendently in each collision cell [32] . 

The strength of the applied pressure field is character- 
ized by the Peclet number, Pe = 'JTi^,, which is here equal 
to the Weissenberg number, where 7 = ggR/{2T]) is the 
shear rate at the cylinder wall and r^j the polymer re- 
laxation time. The Reynolds number Re = gRv„i/rj = 
gR^j / {2r]), where Vm is the maximum fluid velocity, de- 
pends linearly on the shear rate. For the above MPC 
parameters, the viscosity [5S] is such that Re < 1 for all 
considered 7. Equilibrium simulations for a system with 
periodic boundary conditions yield the end-to-end vector 
relaxation time t/j w 3200r. This value agrees within ap- 
proximately 20% with the relaxation time obtained the- 
oretically for a semiflexible polymer with the same ratio 
Lr/Lp [5T]. The relaxation time of the Brownian MPC 
simulation is w 8200t. 

Radial concentration distributions. The radial 
monomer density Pmo is displayed in fig. [5] It is normal- 
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Fig. 2: Radial monomer distributions Pmo for the Peclet num- 
bers Pe = (o), 7 (+), 15 (□), 40 (o), 70 (A), 150 (x), and 
360 (★). Representative error bars are included. Inset: Widths 
of the distributions as function of the Peclet number. 




Fig. 3: Radial center-of-mass distribution functions Pen for 
the Peclet numbers Pe = (o), 7 (+), 15 (□), 40 (o), 70 (A), 
150 (x), and 360 (★). Inset: Widths of the distributions as 
function of the Peclet number with (■) and without (♦) HI. 



ized such that 

/ rP™o(r/a) dr/a^ = 1. (2) 
Jo 

Without flow, there is a depletion zone close to the chan- 
nel wall, which extends approximately one radius of gy- 
ration into the channel. With increasing Peclet number, 
the density in the center of the channel decrease, whereas 
the density near the wall is unchanged. For even larger 
Pe, we observe a migration of the molecule away from 
the surface. At the same time, the density in the channel 
center decreases further and a maximum develops at a fi- 
nite distance from the center. The maximum is shifted to 
smaller radial distances with increasing Pe. Two effects 
contribute to the formation of the maximum. On the one 
hand there is cross-streamline migration of the semiflexible 
polymer due to hydrodynamic interactions. On the other 
hand, the flow field causes an alignment of the molecule 
(see discussion below), which increases the local density 
at radial positions of strong flow-induced alignment. The 
width of the distribution function 

l-R/a 

(r") = / r^Pmo{r/a) dr/a\ (3) 
Jo 

shown in the inset of fig. [U decreases at Peclet numbers 
Pe > 50, refiecting the inward migration of the polymer. 

Figure [3] displays radial center-of-mass distribution 
functions Pcmircm) for various Peclet numbers. Pcm dis- 
plays the same qualitative features as the monomer distri- 
bution Pmo, reflecting the same physical mechanism. For 
Pe < 100, the polymer exhibits a migration away from the 
channel center, which is evident from the decrease of Pcm 
in the central part of the channel, in qualitative agreement 
with experiments |17j . and its increase near the channel 
wall. This is supported by the broadening of the width of 
the distribution with increasing Pe shown in the inset of 
fig. [21 For Pe > 100, the width of the distribution starts 



to decrease with increasing Pe, indicating the presence of 
a wall-induced lift force. As a consequence, the center- 
of-mass distribution shows a clear off-centered peak for 
Pe > 100. 

It is interesting to compare these results with those of 
simulations without HI. We observe an off-center peak in 
the monomer density distribution (compare fig. [5]) also 
without HI, but the difference between the density in cen- 
ter and at the peak position is smaller. However, the width 
(r^) of Pmo does not decrease with increasing Pe in the 
absence of HI. In contrast to the monomer distribution 
functions, the distinct off-center peak of the center-of-mass 
distribution (fig. [3]) is not present in systems without HI. 
Here, the width of Pcm increases monotonically with the 
Peclet number. Hence, without HI, there is an enhanced 
outward migration due to the lack of the wall-induced lift 
force. At a first glance, migration without HI seems sur- 
prising. However, this is related to the suppression of 
steric polymer-wall interactions with increasing Pe due 
to flow- induced polymer alignment (cf. following section). 
Hydrodynamic interactions enhance the migration effect, 
as is seen in fig. [31 Thus, two effects contribute to outward 
migration: alignment of the semiflexible polymer and in- 
tramolecular HI. For the semiflexible polymer, the latter 
dominates. As discussed in ref. [S], HI are less relevant for 
strongly confined fiexible polymers. 

The comparison with confined flexible chains [5] re- 
veals a stronger inward migration for flexible polymers of 
comparable radius of gyration (and thus larger contour 
length), as has already been pointed out in ref. [TB]. The 
reason is the stretching of the flexible polymer along the 
flow direction, which results in a significant longer object 
than the rodlike molecule and hence a more pronounced 
migration. Our studies on the stiffness dependence of mi- 
gration for polymers of equal length show a stronger mi- 
gration for stiffer molecules. 
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Fig. 4: Orientational order parameters S as a function of the 
radial polymer center-of-mass position rem for the Peclet num- 
bers Pe = (o), 7 (+), 15 (□), 40 (o), 70 (A), and 360 
(★). Inset: Probability distribution P of the orientational order 
S = {Scos'^e - l)/2 for the Peclet number Pe = 360 and the 
radial distance rem = 0.5a (■) and rem = 4.5a (♦) 
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Fig. 5: Scaled polymer bending energies Ub as a function of 
radial center-of-mass position for the Peclet numbers Pe = 
(•), 15 (■), 40 (♦), 70 (A), and 360 (★). The increase in 
bending energy near the center is an indication of U-shaped 
conformations. 



Polymer conformations and alignment. In mi- 

crochannel flow with no-slip boundary conditions, the 
shear rate varies linearly across the channel. Hence, the 
force experienced by a polymer depends upon its radial po- 
sition, and its conformations and alignment are expected 
to vary as a function of the radial center-of-mass position. 
Studies of flexible polymers [S] reveal large orientational 
changes by the imposed flow, which is important for cross- 
streamline migration [13j . 

Figure 0] displays the orientational order parameter 

Sir,m.,0)^^{3cosH-l), (4) 

where 6 is the angle between the polymer end-to-end vec- 
tor and the flow direction, as a function of its center-of- 
mass radial position. For Pe < 10, the polymer orienta- 
tion is isotropic in the channel center. For large distances 
from the center, confinement causes a preferred orienta- 
tion along the channel axis; the polymer is almost paral- 
lel aligned near the channel wall. At Pe > 50, flow in- 
duces a strong alignment even at intermediate rem- Most 
importantly, at small radial distances, S decreases again 
for Pe > 100, and becomes even negative for very large 
Peclet numbers, which we attribute to the presence of 
U-shaped conformations (compare fig. [TJi) - as observed 
experimentally [17j . The probability distributions P{S) 
for Tern ~ 0.5a and rem ~ 4.5a reveal the differences in 
alignment. At small rem, the probabilities of U-shapcd 
{S = —0.5) and aligned straight {S — 1) conformations 
are much larger than those of all other orientations. In 
contrast, for larger rem, the probability is largest for al- 
most parallel aligned polymers. We like to emphasize that 
the polymer without HI exhibits very similar orientational 
order parameters, but only at significantly larger Peclet 
numbers. 

The large probability P{S = -0.5) (fig. S]) indicates 



that such conformations are rather stable in the channel 
center. The snapshot of fig.[T|D shows that despite the high 
bending stiffness the polymer adopts U-shapcd conforma- 
tions. Such conformations can be quantitatively analyzed 
by calculating the average bending energy Ub- Figure [5] 
displays Ub as function of the radial center-of-mass po- 
sition and various Peclet numbers. For Pe < 100 the 
bending energy is nearly uniform across the channel cross 
section and decreases near the wall, due to wall induced 
alignment. Its magnitude is close to the thermal aver- 
age Ub = {Lr/a — l)kBT of the nearly harmonic bending 
potential, as expected. An increase in Pe results in an 
increase in its absolute value. For Pe > 100, Ub is dis- 
tinctly higher at the channel center due to the formation 
of a large number of bend conformations. For polymer 
center-of-mass positions slightly out of the center, Ub de- 
creases, since the proportion of bend conformations is re- 
duced. With increasing radial center-of-mass position, the 
local shear rate increases and a polymer assumes transient 
hairpin conformations, as shown in figs.[TJ;,d. The average 
over individual configurations provides a high value for Ub 
for such rem- Adjacent to the channel wall the bending 
energy decreases again. Here, the wall interactions lead to 
a strong alignment with the fiow without any bending. 

Lateral displacement. — The most striking impact 
of flow on the semiflexible polymer is the strong depen- 
dence of its orientation on Pe and the radial distance. This 
aspect provides the key to understand the appearance of 
the off-center maximum in the center-of-mass distribution 
function. 

As is well known, a rod in solution exhibits a larger 
diffusion coefficient parallel to its axis than perpendicular 
to it. In the absence of flow, the semiflexible polymer 
behaves very similar to a stiff rod, we therefore expect 
that the observed polymer orientational differences across 
the channel will lead to differences in the lateral dynamic 
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Fig. 6: Mean radial polymer center-of-mass displace- 
ments (Ar) = ([rcm(to + At) - rcm(io)]ecm(io)> as func- 
tion of its radial center-of-mass position rem (to) for the 
various indicated Peclet numbers. Inset: Radial poly- 
mer center-of-mass mean square displacements (Ar^) — 
{{rcm{to + At)ecm(to) - {r^mito + At)ecm(to)>)^)- 



translational behavior. 

In order to characterize this behavior, we present in 
fig. ini the radial center-of-mass displacement (Ar) = 

(Km(io+At)-rc,„(io)]ecm(io))), Gem = T"cm/kcm|, with 

respect to the initial radial position rcm.(io) of the polymer 
within a time At, which is chosen such that the center-of- 
mass displacement is less than a monomer size a and the 
mean square displacement is linear in time. Figure [S] shows 
a drift away from the wall for all Peclet numbers and an 
outward drift for the largest Pe in the channel center, in 
agreement with the center-of-mass distribution functions 
of fig. |3l At Pe = 0, the drift is governed by steric rod- 
wall interactions for rem ^ 1-5, which is reflected by the 
decrease of Ar for that range of Vcm ■ An increase of the 
flow rate reduces the steric effect by aligning the polymer, 
however, the wall lift force now leads to an increase of the 
magnitude of the drift for rem/ a ^ 5. 

The mean square displacement ^Ar^^ 

{{rcm{ta + At)ecm(io) - {Vcmih + At)ecm{to)))^), 
shown in the inset of fig. [51 displays only a slight 
dependence on drift for rc^n/a ^ 5; its decrease with 
increasing rem refiects the dependence of the mean 
square displacement on radial position. Without fiow, 
the polymer orientation is isotropic in the central part of 
the channel, and it is essentially parallel to the channel 
axis near the wall. Hence, we see the difference between 
isotropic and perpendicular diffusive motion originating 
from intramolecular hydrodynamic interactions. For 
larger radii, steric wall interactions lead to a decrease of 
the mean square displacement for low Peclet numbers, 
whereas lift forces increase (Ar^) for larger Pe values. 

Non-hydrodynamic simulations yield constant mean 
square displacements for rem /a ^ 5 for the considered 
range of Peclet numbers. At larger r^mi steric wall inter- 
actions lead to a decrease of (Ar^), similar to that for 



Pe = in fig. H 

Hence, we conclude that the center-of-mass distribution 
function is determined by a competition of the radial out- 
ward diffusion of a flow-aligned rod and inward migration 
caused by a wall-induced lift force [T^ . 

Polymer tumbling. — Flexible polymers in shear 
and Poiseuille flow undergo large conformational changes 
by tumbling motion, i.e., the polymers exhibit a cyclic 
stretching and collapse dynamics [2T1[33]. Visual exam- 
ination of the dynamics of the semiflexible polymer also 
indicates tumbling motion. This cyclic motion is charac- 
terized by the tumbling time It- Analytical calculations 
for a planar shear flow have shown that tr ^ 'y~2/3 
flexible and rodlike polymers, where 7 is the shear rate 
|21| . In order to compute the tumbling time for different 
Pe, we determine the cross-correlation function between 
the mean extensions 5z = z — (z) and Sr = r — (r) in flow 
and radial directions, respectively, which is defined as 

C.rit) = {6zito)Srito + t))/.J{6z^ito)}{Sr'{to)). (5) 

The time dependence of the correlation function is pre- 
sented in fig. [T] MicroChannel flows are more complex in 
nature compared to planar shear flows. In planar shear 
flow, the polymer exhibits essentially a two-dimensional 
behavior, whereas the channel flow is rotational symmet- 
ric with respect to the channel axis and more complex 
conformational transitions appear. Therefore, Czr is qual- 
itatively different from to a similar correlation function in 
planar shear flow for flexible polymers [S^- In the latter 
case, at Pe > 1, the correlation function is asymmetric 
and displays a minimum and a maximum, whereas we ob- 
tain a time-symmetric negative correlation function with 
a peak at t = 0. The width of the correlation function de- 
creases with increasing flow strength, which is a measure 
for the dependence of tx on Pe. Hence, we quantify the 
tumbling time by the full width at half maximum of C'zr 
for each value of Pe. As shown in the inset of fig. [71 the 



_ tumbling time decreases as tx ~ Pe 



-0.66 



i.e., the expo- 



nent is very close to the analytically predicted exponent 
for planar shear flow [21]. In order to verify the result, a 
tumbling time is also extracted from the distribution func- 
tion of the time intervals between consecutive crossings of 
the end-to-end vector with the plane perpendicular to the 
flow direction [5]. This distribution decays exponentially 
{cxp{—t / t'rp)) for sufficiently large times. The tumbling 
time t'rp exhibits the same dependence on the Peclet num- 
ber as tx obtained from C'zr] quantitatively it is smaller 
by approximately a factor 1.5. Hence, both definitions of 
tumbling time are equivalent for the considered system. 

Summary and Conclusions. We have analyzed 
the flow behavior of an actin-like semiflexible polymer con- 
fined in a cylindrical channel. The Poiseuille flow induces 
a strong radial-dependent polymer alignment parallel to 
the channel associated with a radially outward migration 
of its center of mass at low Peclet numbers [T31[TS] and 
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Fig. 7: Cross-correlation functions Czr for the Peclet num- 
bers Pe = 15, 40, 70, 150, 360, and 570 (bottom to top). Inset: 
Tumbling times tx calculated from the cross-correlation func- 
tion Czr (B) and the probability distribution of time intervals 
between radial orientations of the end-to-end vector (♦). 



an inward migration at large Pe; both effects are due to 
hydrodynamic interactions. 

This behavior is surprisingly similar to that of a flex- 
ible polymer confined in a channel pi[^[51[51-[TT|. It is 
caused by the common (flow induced) anisotropy of the 
polymers, which is responsible for the appearance of mi- 
gration jl3j . Flexible polymers, however, exhibit much 
more pronounced conformational changes and large defor- 
mations, which depend on the radial center-of-mass posi- 
tion, but there is no enhanced probability for U-shaped 
conformations. We also find different dependencies of the 
tumbling time on the flow strength [^j • 

Our studies confirm the importance of HI in microchan- 
nel flows. Nevertheless, astonishingly similar distributions 
and conformations are obtained in simulations without 
HI, but they appear at significantly larger Peclet num- 
bers. However, in this large Pe regime, hydrodynamics 
determines the distributions in the channel. Thus, non- 
hydrodynamic simulations may provide reasonable struc- 
tural and conformational properties of a fully hydrody- 
namic system as long as wall lift forces are negligible, but 
will fail for large Pe even qualitatively. 
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